The Antiferromagnetic Heisenberg Model on Clusters with 

Icosahedral Symmetry 

N.P. Konstantinidis* 
Department of Physics and Department of Mathematics, 
University of Dublin, Trinity College, Dublin 2, Ireland 
(Dated: February 2, 2008) 

Abstract 

The antiferromagnetic Heisenberg model is considered for spins s, = i located on the vertices of 
the dodecahedron and the icosahedron, which belong to the point symmetry group J^. Taking into 
account the permutational and spin inversion symmetries of the Hamiltonian results in a drastic 
reduction of the dimensionality of the problem, leading to full diagonalization for both clusters. 
There is a strong signature of the frustration present in the systems in the low energy spectrum, 
where the first excited states are singlets. Frustration also results in a doubly-peaked specific 
heat as a function of temperature for the dodecahedron. Furthermore, there is a discontinuity 
in the magnetization as a function of magnetic field for the dodecahedron, where a specific total 
spin sector never becomes the ground state in a field. This discontinuity is accompanied by a 
magnetization plateau. The calculation is also extended for Si = 1 where both systems again 
have singlet excitations. The magnetization of the dodecahedron has now two discontinuities in an 
external field and also magnetization plateaux, and the specific heat of the icosahedron a two-peak 
structure as a function of temperature. The similarities between the two systems suggest that 
the antiferromagnetic Heisenberg model on a larger cluster with the same symmetry, the 60-site 
cluster, will have similar properties. 

PACS numbers: PACS numbers: 75.10.Jm Quantized Spin Models, 75.50.Ee Antiferromagnetics, 75.50.Xx 
Molecular Magnets 
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I. INTRODUCTION 



The antiferromagnetic Heisenberg model has been the object of intense investigation 
in recent years as a prototype of strongly correlated electronic behavior 0- The effects 
of low dimensionality, quantum fluctuations and frustration combine together to produce 
new phases different from conventional Neel-like order, where the order parameter is not the 
staggered magnetization or there is a lack of a local order parameter altogether 0,0,0). They 
can also have dramatic consequences on the energy spectrum, as is the case for frustrated 
multiple spin exchange models and the Kagome lattice, where the low energy excitations 
are singlets [5j. Specific heat calculations of small Kagome lattice samples have revealed a 
two-peak structure, with the first peak below the singlet-triplet gap Jo, 13 - The double-peak 
structure was also obtained for other multi-spin exchange models a |9( , a pyrochlore slab 
[loj | . and the A chain [llj]. It is a natural question to ask whether there are other frustrated 
systems with such unconventional behavior. 

Here clusters with the connectivity of the fullerenes will be considered QQ. These 
molecules weakly bound to form crystals which become superconductors when doped with 
alkali metals |l4l . Il5j | . Their transition temperature is above 40 K, a transition tempera- 
ture much higher than the one of conventional superconductors. Chakravarty and Kivelson 
suggested that an electronic mechanism at intermediate scales is responsible for supercon- 
ductivity in Ceo, the fullerene with 60 carbon atoms, when doped with alkali metals [16] . 
It is an open question if a repulsive interaction can produce pairing in such systems. The 
Hubbard model has been used to investigate an electronic mechanism for superconductiv- 
ity on this cluster 0, An exact treatment of the model in the full Hilbert space of 
the molecule is prohibitive due to its size. Therefore, as a first step, smaller molecules of 
the fullerene type were considered to gain insight, as well as the strong on-site repulsion 



limit of the Hubbard model at half-filling, the antiferromagnetic Heisenberg model 
Coffey and Trugman found that connectivity and frustration lead to non-trivial behavior 
at the classical level in a magnetic field To study the effect of quantum fluctuations 
on the classical results a 20-site cluster, the dodecahedron, was considered (figure [T]). It 
is three-fold coordinated, has all sites equivalent and consists of 12 pentagons. The model 
was studied with perturbation theory around the classical limit and was found to possess 
a singlet as the first excited state, and a discontinuity in the magnetization as a function 
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of quantum fluctuations [2l|. This discontinuity was originally found at the classical level 
|l9| . Therefore it is of interest to study the structure of the low energy spectrum of the 
dodecahedron and its response in a magnetic field for Sj = \. A more general question is if 
the presence of singlets in the excitation spectrum is also a property of other clusters, and 
if there is a correlation of magnetic behavior with space group symmetry and connectivity 



22| . The determination of the full energy spectrum will reveal if the double-peak specific 



heat structure is also a property of the dodecahedron. However, the projection of the total 
spin on the z axis, S z , is not a good quantum number when perturbing around the classical 
state. In addition, multiprecision arithmetic was needed to analytically continue the series 
expansions in j^jj. Thus the full calculation of the energy spectrum was prohibited due to 
memory requirements, even though symmetry was partially taken into account. The ground 
state properties of the model were firstly calculated by Modine and Kaxiras 2^. 

In this paper the antiferromagnetic Heisenberg model is studied on the dodecahedron and 
a smaller cluster with 12 equivalent sites and the same spatial symmetry, the icosahedron 
figure , for Sj = | and 1. Permutational and spin symmetries are taken into account 

ization except for the dodecahedron when Sj = 1, 
28|. The symmetry group of the clusters is 1^, the 



24 |2jj |2y], and this leads to full diagona 



where Lanczos diagonalization is used 
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largest point symmetry group with 120 operations |29j. The icosahedron has 20 triangles, 
and the spins on the vertices are five-fold coordinated. Therefore the two clusters share the 
same symmetry but their connectivity is different. The results for the low energy spectrum 
of the dodecahedron for s» = | show that the ground state is a singlet, while the first excited 
state is also a singlet but five-fold degenerate. The next excited state is also a singlet, and 
then a series of states with total spin S = 1 follows. This series of low-lying singlets is a 
consequence of the connectivity and frustration of the model. The specific heat also shows a 
non-conventional behavior as a function of temperature, with a peak inside the energy gap 
and a second peak at higher temperatures. The low energy spectrum of the icosahedron is 
similar with the one of the dodecahedron, with the same ordering and relative spacing of 
low energy levels. The specific heat has a well-pronounced peak inside the energy gap, but 
there is no second peak but rather a shoulder at higher temperatures. 

For spins with magnitude Sj = 1, the ground and first excited state of the dodecahedron 
are closely spaced singlets. The next excited state is a three-fold degenerate triplet. The low 
energy spectrum for the icosahedron is again the same in the symmetry, the ordering and 
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the relative spacing of the low lying levels. The specific heat can only be calculated for the 
icosahedron and has now two peaks, however there is no peak inside the energy gap. This 
along with the reduced number of low energy singlet states compared with the full quantum 
Si — | case, indicates a change in the spectrum with increasing Sj. 

The behavior of the magnetization in a magnetic field is non-trivial for the dodecahedron. 
There is a discontinuity as the energy of a particular total spin S sector never becomes the 
ground state in a field. This feature, observed at the classical level by Coffey and Trugman 



19j . survives in both the = ~ and 1 cases, twice in the latter. Similar behavior is not 



observed for the icosahedron, which also has discontinuous magnetization in a field at the 
classical level. Even though the two clusters have the same spatial symmetry, the behavior 
in a field appears to depend on their polygon structure. The magnetization discontinuities 
in the dodecahedron are accompanied by magnetization plateaux. 

The similarities of the dodecahedron and icosahedron spectra suggest that predictions 
about the low energy structure can be made for a larger cluster of fullerene type connectivity 
with the same symmetry J/j, the 60-site cluster, where again all sites are equivalent Even 
though the 12- and 20-site clusters have different coordination number and consist of different 
types of polygons, the structure of the low energy spectrum is the same. The 60-site cluster 
has also a discontinuity in the magnetization in a field at the classical level [l9|. However, 
it is not clear what the response in a magnetic field will be for the quantum case, since the 
two smaller clusters have different behavior, which seems to depend on the connectivity as 
well as their polygon structure. Similar considerations could correlate the behavior of more 
complicated models with orbital degrees of freedom like the Hubbard model on the 60-site 
cluster with the behavior of the same models on the dodecahedron and the icosahedron. 

The plan of this paper is as follows: in section [H] the model and method are introduced, 
in section IIHI the low energy spectra of the two clusters are presented for both Sj = ~ and 
1, in section HYI specific heat and magnetic susceptibility data are presented, in section El 
the ground state magnetization is considered and finally in section IVII the conclusions are 
presented. 
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II. MODEL AND METHOD 



The Hamiltonian for the antiferromagnetic Heisenberg model is 

H = J si • Sj - hS z (1) 

<i,j> 

where the spins are located at the vertices of the clusters and < i,j > denotes nearest 
neighbors. The coupling constant J is positive and will be set to 1 from now on, defining 
the unit of energy, h is the strength of an external magnetic field. 

Minimization of memory requirements for diagonalization is possible with the use of the 
symmetries of the model. These include permutational and spin space symmetries [24. I25I 
l26j |. The Hamiltonian commutes with S and S z . However, even though it is straightforward 
to work in an S z subspace, there is no efficient method to construct symmetry adapted 
eigenstates of S. The Hamiltonian is symmetric under combinations of permutations of the 
spins that respect the connectivity of the cluster. The group of permutations is the symmetry 
group of the cluster in real space j^J. The model also posseses time- reversal symmetry and 
inverting the spins is a symmetry operation in the S z = sector. The corresponding 
group is comprised of the identity and the spin inversion operation. The full symmetry 
group of the Hamiltonian is the product of the space group and the group of spin inversion. 
Taking the full symmetry into account, the S z basis states can be projected into states that 
transform under specific irreducible representations of the symmetry group. In this way the 
Hamiltonian is block-diagonalized into smaller matrices, and their maximal dimension is 
dramatically reduced compared to the full Hilbert space size. The largest sub-matrix of the 
block-diagonalized Hamiltonian for the dodecahedron has dimension 7, 058 for spins s« = |, 
therefore full diagonalization is possible. For s» = 1 only a few of the lowest eig envalues 
for each irreducible representation were obtained with Lanczos diagonalization j^j , and 
for the subspace S z = 1 this was not possible for the five-dimensional representations. The 
largest matrix for which the lowest eigenvalue was found has dimensionality 13, 611, 598 and 
is complex. In the case of the icosahedron full diagonalization is possible for both Sj = | 
and 1, the largest matrix having a dimension of 2, 982 in the latter case. 
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III. LOW ENERGY SPECTRA AND CORRELATION FUNCTIONS 



The low energy spectrum of the icosahedron is presented in table The ground state 
energy is a singlet and belongs to the irreducible representation A u>s , where the first index 
g or u denotes symmetry or antisymmetry with respect to space inversion 2^|, and the 
second s or a symmetry or antisymmetry with respect to spin inversion when S z = 0. The 
energy per spin equals —0.51566. The first excited state is also a singlet, belongs to the 
representation H g s and is five-fold degenerate. The next excited state is also an S = 
state, it is non-degenerate and belongs to A 9j8 . Then a series of S = 1 states follows. The 
gap to the first excited state is 0.53344 and the singlet-triplet gap 0.89988. 

The nearest neighbor correlation function for the ground state equals —0.20626, almost 
four times less in strength than the value of an isolated dimer, —0.75. For the first excited 
state there are two different types of nearest neighbor correlation functions. The first equals 
-0.16930 for the pairs (1, 3), (9, 11), (5, 8), (6, 7), (2, 10), (4, 12) in figure 121 and the second 
—0.19328 for the rest of the bonds. The six pairs face each other in pairs and belong in 
different triangles. There are five different ways of distributing the pairs like that on the 
icosahedron, thus the state is five-fold degenerate. The second excited singlet has all nearest 
neighbor correlation functions equal to —0.18748. The correlation functions other than the 
nearest neighbor in the ground state between site 1 and the rest of the sites are equal to 
—0.13908 for spin 11, more than half the value of the nearest neighbor correlation function, 
and 0.08408 for the rest of the spins. 

The low energy spectrum of the dodecahedron is shown in table |H] and has the same 
structure as the one of the icosahedron, after comparison with table |U In particular, the 
four lowest energy states belong to exactly the same irreducible representations as the cor- 
responding states in the icosahedron. The spacing of the lowest energy states is also similar. 
The ground state energy equals —9.72219 |2ll.l23j|. and the energy per spin in the ground state 
is now higher and equal to —0.48611. Classically, the ground state energy per bond is — ^ 
for the dodecahedron |2]J and — ^ for the icosahedron (previously reported by Schmidt and 
Luban in j^), therefore the energies per spin are equal. Quantum fluctuations reduce the 
energy more for the cluster with the highest coordination number, even though the number 
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of bonds is the same for both. The gap to the first excited state is also smaller and equal to 
0.31567, and the same is true for the singlet-triplet gap which is 0.51383. The nearest neigh- 
bor correlation function is equal to —0.32407 for the ground state quite stronger in 
magnitude than the one of the icosahedron, as in the classical case. This is attributed to 
the lower coordination number of the dodecahedron. For the first excited state, the value 
is -0.33585 for the pairs of spins numbered (1,2), (9, 10), (14, 15), (4, 12), (17, 18), (7,20) in 
figure ^ while for the rest of the pairs it is —0.30797. The above six pairs are facing each 
other in pairs in the dodecahedron and belong in different pentagons. As in the icosahedron, 
there are five different ways of distributing the six pairs of spins in the above manner, thus 
the state is five-fold degenerate. The bonds on these pairs are more singlet like than the 
ground state, unlike the icosahedron. Another difference is that the bonds on the six pairs 
have now lower energy than the rest of the bonds. For the next excited state, the nearest 
neighbor correlation functions are equal to —0.31175. The correlation functions other than 
the nearest neighbor in the ground state between site 1 and the rest of the sites are equal 
to 0.06540 for spins 3, 4, 7, 8, 14 and 15, -0.03882 for spins 9, 10, 12, 13, 16 and 20, 0.03307 
for spins 11, 17 and 19 and —0.03649 for spin 18. They are significantly smaller than the 
nearest neighbor correlations. This is in contrast to the icosahedron, where correlations 
survive longer distances. 

B. Si = 1 

Full diagonalization is still possible for the icosahedron when Sj = 1. The low energy 
spectrum is shown in table IIHI The ground state is a singlet in the A g>s irreducible repre- 
sentation and its energy is —18.56111. It is now symmetric with respect to space inversion, 
in contrast to the Si = ~ case. The first excited state is a singlet in the A u>s representation, 
which included the ground state for = |, with energy close to the ground state and equal 
to —18.42539. The next excited state is a three-fold degenerate triplet in the representation 
T 2uA with energy —17.83998, followed by another three- fold degenerate triplet in the repre- 
sentation T 2 g t a with energy —17.80499. The energy per spin in the ground state is —1.54676 
and the nearest neighbor correlation function —0.61870. This value is to be compared with 
the singlet ground state energy of a dimer with spins Sj = 1, which equals —2. The ratio 
of the two values is larger compared with the corresponding ratio for the Sj = ~ case and 
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closer to the classical result. For the first excited state the nearest neighbor correlation is 
—0.61418, much closer to the ground state value compared with the Sj = | case. The gap 
to the first excited state equals 0.13572 and the singlet-triplet gap 0.72113. The correlation 
functions other than the nearest neighbor in the ground state between site 1 and the rest of 
the sites are equal to —0.74630 for spin 11, and 0.36796 for the rest of the spins. Compared 
with the Si = \ case, the magnitudes of the next than nearest neighbor correlation functions 
are significantly increased with respect to nearest neighbors, and the correlation with site 
11 is even stronger than the nearest neighbor correlation. 

As was the case for = |, the low energy spectrum (table HV]) is similar for the clusters 
when Si = 1. The ground state for the dodecahedron is a non-degenerate singlet with 
energy —30.24551, and an energy per spin equal to —1.51228. Again the energy per spin is 
higher than the one of the icosahedron. The first excited state is also a singlet close to the 
ground state with energy —30.21750, and the next two excited states are triply degenerate 
triplets with very close energies, —29.92161 and —29.91011. The energy gap is 0.02801 
and the singlet-triplet gap 0.32390, values again smaller than the ones of the icosahedron. 
The nearest neighbor correlation functions in the ground and first excited states are equal 
to —1.00818 and —1.00725 respectively. The correlation functions other than the nearest 
neighbor in the ground state between site 1 and the rest of the sites are equal to 0.27896 for 
spins 3, 4, 7, 8, 14 and 15, -0.20912 for spins 9, 10, 12, 13, 16 and 20, 0.35962 for spins 11, 17 
and 19 and —0.47333 for spin 18. Similarly to the icosahedron, the magnitudes are increased 
with respect to the nearest neighbor correlations compared with the s« = \ case. Farther 
than nearest neighbor correlations are stronger for the icosahedron, as was the case for 
Sj = ~. It wasn't possible to diagonalize the five-fold degenerate irreducible representations 
for S z — 1. Therefore the 5* values in parentheses in table HVl are deduced by comparison 
with the states of the icosahedron, since the spectra are similar. In any case these S values 
can only be or 1, with these energies absent from the S z = 2 spectrum. 

IV. SPECIFIC HEAT AND MAGNETIC SUSCEPTIBILITY 

The temperature dependence of the specific heat and the magnetic susceptibility for the 
cases where exact diagonalization is possible is shown in figures El and |U respectively. For 
the icosahedron and s, = ~, there is a peak in the specific heat around T = 0.219 and a 
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shoulder around T = 0.8. The peak is inside the energy gap, a feature characteristic of 
frustrated systems 3, El- For the dodecahedron there are two well defined peaks. The first 
peak is centered around T = 0.120 and the second around T = 0.627. The first peak is 
inside the energy gap. Similar results have been obtained for other frustrated systems, and 
Sindzingre et al. attributed the peak to the combined effect of singlet excitations inside the 
singlet-triplet gap and low-lying triplet excitations for the Kagome lattice p|. In contrast, 
Syromyatnikov and Maleyev considered a Kagome star where the peak results from the 
increase in the density of states just above the spin gap In the present case both 
the 5 = and 5=1 sectors contribute for the peak inside the energy gap. The 5 = 
contribution is due to the low-lying singlets, while the 5 = 1 contribution is also due to a few 
of the lowest energy triplets. For the shoulder and the higher energy peak in the two systems 
respectively higher 5 sectors contribute as well. For the icosahedron and Sj = 1, there are 
three regions in the graph. At temperatures lower than the energy gap, the specific heat 
initially rises slowly towards an area between T pa 0.055 to 0.07 where it almost stabilizes 
with temperature. The contributions come from the lowest energy states in the 5 = and 
5 = 1 sectors. The specific heat then increases rapidly reaching a peak around T = 0.395, 
followed by a slight decrease to a local minimum at a temperature around the singlet-triplet 
gap, T = 0.619. Following that it increases again to another peak around T = 0.786 and then 
starts dropping towards zero. The contributions to the two peaks and the local minimum 
come from various spin sectors. Even though the double-peak structure is similar to the 
specific heat of the dodecahedron for s« = |, there is no peak inside the energy gap, and this 
signifies a change with respect to the Sj = \ case. The magnetic susceptibility is plotted 
in figure HI with a temperature dependence similar for the three cases. The position of the 
maximum follows the pattern for the lowest temperature peak of the specific heat, with a 
peak temperature decreasing with size and increasing with spin magnitude. 

V. GROUND STATE MAGNETIZATION 

Frustrated spin systems have been found to exhibit magnetization jumps in the presence 
of an external field, where the lowest state in a particular 5 sector never becomes the ground 
state in the field [3l|. The magnetization in an external field has been calculated for the 
dodecahedron at the classical level, and such a discontinuity was found A discontinuity 
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in the magnetization was also found when the classical state is perturbed with the quantum 

ri 

fluctuations [21]. The calculation for the classical ground state of the icosahedron shows also 
the presence of a magnetization jump in a field. In the quantum case, the lowest energy in a 
magnetic field in a particular S sector is calculated from the lowest energy in the absence of a 
field by adding the Zeeman term. The lowest energies in the different S sectors are shown in 
tables IVl and IVII for the dodecahedron for Sj = | and s« = 1 respectively. As was mentioned 
before, calculation of the lowest energies is not possible for the five-dimensional irreducible 
representations when S z = 1. However, by comparing the spectra of the five- dimensional 
representations for S z = with the spectra of the rest of the representations for S z — 1, it is 
seen that the lowest energy for S = 1 does not belong to the five-dimensional representations 
(table HVj) . The ground state magnetization is found by comparing the energies in a field in 
the different sectors. The plots for the dodecahedron for both Sj = | and 1 are shown in 
figure |3J where the magnetization M is the total spin S normalized to the number of sites 
and the magnitude of each spin, with steps between different S sectors equal to 0.1 for Sj = | 
and 0.05 for Sj = 1. There is a discontinuity of M between 0.4 and 0.6 for Sj = ~ (S = 4 and 
6), and two discontinuities for s« = 1, between 0.4 and 0.5 (S = 8 and 10), and between 0.7 
and 0.8 (S = 14 and 16). These discontinuities are associated with magnetization plateaux 
on both sides of the jumps j^. Magnetization plateaux are also observed in figure El for 
M = when s< = |, and for M = 0.2 and 0.3 for s t = 1. 

As seen from tables and IVI[ the ground states on each side of the jumps are non- 
degenerate, and the spatial symmetry of the ground state switches from symmetric to anti- 
symmetric as the field increases. Spin inversion symmetry has only be determined for s» = \ 
for the cases of interest, and in that case this symmetry does not change as the discontinuity 
takes place. The lowest state of M that never becomes the ground state is in every case 
symmetric with respect to spatial inversion, and it belongs to the same irreducible repre- 
sentation for the Si — ^ an d the lowest field Sj = 1 jump. It is also interesting that in the 
Si — 1 case, even though the lowest state for S = 6 is non-degenerate as is the S — 8 state, 
there is no discontinuity associated with these two states, indicating the significance of the 
different spatial inversion symmetry for the states on the two sides of a magnetization jump. 
However, a magnetization plateau exists for the corresponding magnetization M = 0.3. 

The behavior of the correlation functions on either side of the magnetization jumps is 
shown in table IV11I The nearest-neighbor and next nearest-neighbor correlation functions, 
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as well as Si ■ Sn, become more positive as the magnetic field is getting stronger. For Si • Sis, 
the magnetization jump decreases its value with increasing field. However, this weakening 
is very pronounced for the S = 8 to S = 10 transition for Sj = 1. For Si ■ Sg, the lower 
jump in the Sj = 1 case decreases its strength, contrary to what happens for the other two 
discontinuities. It is concluded that the first discontinuity in the Sj = 1 case has different 
characteristics from the other two jumps as far as longer range correlation functions are 
concerned. This could possibly be due to the different change in behavior with respect to 
spin inversion as this discontinuity takes place, compared with the other two discontinuities. 

The data for the two lowest quantum numbers Sj shows that the discontinuity in the 
classical solution survives the quantum fluctuations. The presence of jumps for s» = |,1 
and oo indicates that this is probably a characteristic of the system for any quantum number. 
The common features of the mechanism of the effect for the two lowest quantum numbers 
also point in this direction. For the icosahedron there is no such jump for s,- = ~ or 1, even 
though there is a jump at the classical level, and quantum fluctuations are stronger. It is 
therefore concluded that the discontinuity in the magnetization at the quantum level is not 
a consequence of the symmetry, but rather of the polygon structure of the clusters. For 
larger clusters with non-equivalent sites and including hexagons no such discontinuity was 



found 



22j . Thus the pentagon-only structure of the dodecahedron appears crucial for the 



presence of the discontinuities. 



VI. CONCLUSIONS 



The antiferromagnetic Heisenberg model was solved on the icosahedron and the dodec- 
ahedron, which belong to the point symmetry group Ih- It was found that the low energy 
spectra are similar, and consist of excited states that are singlets for both — | and 1. 
Spin correlations survive longer distances for the icosahedron. The specific heat calculation 
revealed a non-trivial dependence on temperature, with more than one peak for the dodeca- 
hedron for Si = |, and the icosahedron for Sj = 1. This feature has also been found for other 
frustrated spin systems 0]-^^. There are discontinuities in the ground state magnetization 
as a function of magnetic field for the dodecahedron for both Sj = | and 1, but this is not the 
case for the icosahedron. This discontinuity was found at the classical level [3], and here 
it has been shown to survive the quantum fluctuations. It appears to be a characteristic of 
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the system for any quantum number. There is also a number of magnetization plateaux for 
both quantum numbers. These non-trivial properties of the clusters are a consequence of 
the frustrated interactions and the connectivity. Other clusters of fullerene-type geometry 
show similar properties in their low energy spectrum 2^. The similarity in the spectra 
and the specific heat of the icosahedron and the dodecahedron points in similar properties 
for the larger cluster with the same symmetry and 60 sites, also equivalent. It also shows 
that more complicated fermionic models for the 60-site system, such as the Hubbard model, 
could instead be considered on the icosahedron and the dodecahedron, and the results will 
be reliable predictions for the properties of the models on the larger system. However, 
no prediction can be made for the response in a field for the antiferromagnetic Heisenberg 
model, which depends not only on the symmetry but also on the polygons that make up the 
clusters. The dodecahedron exhibits a jump and has only pentagons, while the 60-site sys- 
tem has hexagons as well. It has been found that hexagon correlations are more important 
for frustrated systems of the fullerene type with more than 20 sites, for which there is no 
magnetization discontinuity 22| . Another point deserving further attention is the transition 
from the quantum to the classical limit, regarding the low energy spectra and the presence 
of singlets, the temperature dependence of the specific heat and the presence of magnetiza- 
tion discontinuitites. The general behavior of the clusters with Ih point group symmetry is 
determined to a significant extent by the symmetry, however the coordination number and 
the polygons that make up the structures are also important for the determination of their 
properties. 
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TABLE I: Low energy spectrum for the icosahedron for Sj = i. The first index of the irreducible 
representations indicates the behavior under spatial inversion, where g stands for symmetric and 
u for antisymmetric 



29], and the second under spin inversion for the S z = component, where s 
stands for symmetric and a for antisymmetric. 



energy 


multiplicity 


irreducible 


S 


energy 


multiplicity 


irreducible 


s 






representation 








representation 




-6.18789 


1 







-4.82887 


5 







-5.65445 


5 


H 9,s 





-4.76398 


5 




1 


-5.62426 


1 


Ag,s 





-4.50000 


4 


P 





-5.28801 


3 


T 2g ,a 


1 


-4.50000 


1 







-5.17728 


3 




1 


-4.44972 


4 


7? 

± u,a 


1 


-5.10989 


3 


T 2u ,a 


1 


-4.31691 


5 


H 9,a 


1 


-4.86352 


4 


/? 

r 9,a 


1 


-4.29999 


3 


T2 u ,a 


1 



TABLE II: Low energy spectrum for the dodecahedron for Sj = |. Notation as in table □ 



energy 


multiplicity 


irreducible 


s 


energy 


multiplicity 


irreducible 


S 






representation 








representation 




-9.72219 


1 


A 





-8.87964 


5 







-9.40651 


5 


H g,s 





-8.69499 


3 


T 2u ,a 


1 


-9.35236 


1 







-8.69441 


5 


H u ,a 


1 


-9.20836 


3 


T 2 g,a 


1 


-8.66571 


3 




1 


-9.18649 


4 


7? 

1 u,a 


1 


-8.65030 


4 


7? 


2 


-9.13048 


3 


T 2u ,a 


1 


-8.63614 


5 


H g,a 


1 


-8.97112 


3 


Tlg,a 


1 


-8.63178 


4 


7? 


1 
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TABLE III: Low energy spectrum for the icosahedron for s, = 1. Notation as in table □ 



energy 


multiplicity 


irreducible 


S 


energy 


multiplicity 


irreducible 


s 






representation 








representation 




-18.56111 


1 







-17.15212 


3 




1 


-18.42539 


1 


A 





-17.00416 


5 


Hg,s 





-17.83998 


3 


T 2u ,a 


1 


-16.97485 


4 


p 

± u,a 


1 


-17.80499 


3 


T 2g ,a 


1 


-16.82345 


5 


H g,a 


1 


-17.60137 


5 







-16.75333 


5 




1 


-17.19717 


4 


p 





-16.74705 


1 


Ag,s 





-17.17453 


5 


H g,s 





-16.73750 


4 


P 


1 



TABLE IV: Low energy spectrum for the dodecahedron for Sj = 1. Notation as in table D The S 
numbers in parentheses can be or 1 since they are missing from the S z = 2 spectrum, and they 
are assigned values after comparison with the icosahedron spectrum in table ITTT1 



energy 


multiplicity 


irreducible 


s 


energy 


multiplicity 


irreducible 


S 






representation 








representation 




-30.24551 


1 


Ag,s 





-29.61145 


3 




1 


-30.21750 


1 


A 





-29.57332 


4 


p 





-29.92161 


3 


T 2 u,a 


1 


-29.50512 


3 




1 


-29.91011 


3 


T 2g ,a 


1 


-29.45063 


4 




1 


-29.85881 


5 


H g,s 


(0) 


-29.39457 


5 




(1) 


-29.67223 


5 




(0) 


-29.35464 


5 


H 9,a 


(1) 


-29.65951 


4 


p 

± u,a 


1 


-29.31754 


3 


T 2u ,a 


1 
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TABLE V: Lowest energy in the various S sectors for the dodecahedron for Sj = |. Notation as 
table □ 



s 


energy 


multiplicity 


irreducible 


S 


energy 


multiplicity 


irreducible 








representation 








representation 





-9.72219 


1 




6 


-2.47099 


1 




1 


-9.20836 


3 


T 2 g,a 


7 


-0.13397 


3 




2 


-8.65030 


4 


77 


8 


2.34152 


5 


H 9,s 


3 


-7.72967 


4 


77 

r 9,a 


9 


4.88197 


3 


T2u,a 


4 


-6.40730 


1 


Ag,s 


10 


7.5 


1 




5 


-4.37206 


4 


77 

r 9,a 





TABLE VI: Lowest energy in the various S sectors for the dodecahedron for Si = 1. Notation 
in table [fl For sectors with S higher than 1 the symmetry under spin inversion for the S z = 
component has not been determined (except for S = 20). 



s 


energy 


multiplicity 


irreducible 


S 


energy 


multiplicity 


irreducible 








representation 








representation 





-30.24551 


1 


Ag,s 


11 


-9.72811 


3 




1 


-29.92161 


3 


T 2u ,a 


12 


-6.26689 


5 


H u 


2 


-29.30598 


5 


H 9 


13 


-2.61717 


4 


F 9 


3 


-28.39788 


3 


T lu 


14 


1.22128 


1 


A 9 


4 


-27.20612 


4 


F 9 


15 


5.52324 


3 


T2 9 


5 


-25.49868 


4 


F 9 


16 


9.80295 


1 


A u 


6 


-23.64370 


1 


A 


17 


14.65742 


3 


Tiu 


7 


-21.26717 


4 




18 


19.65432 


5 


H 9 


8 


-18.82343 


1 




19 


24.76393 


3 


T 2u 


9 


-15.88505 


4 




20 


30 


1 


Ag,s 


10 


-12.98367 


1 
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TABLE VII: Unique ground state correlation functions for the lowest energy states on either side 
of the magnetization discontinuities. 



correlation 


„ _ 1 
b i — 2 


„ _ 1 
b i — 2 


Si — 1 


Si — 1 


Si — 1 


Sj — 1 


function 


S = 4 


S = 6 


S = 8 


S = 10 


S = 14 


S = 16 


Si ■ S2 


-0.21358 


-0.08237 


-0.62745 


-0.43279 


0.04071 


0.32677 


Si-S 3 


0.08764 


0.12061 


0.34296 


0.43033 


0.60512 


0.72748 


Si ■ Sg 


0.04971 


0.08038 


0.18218 


0.16926 


0.47185 


0.59169 


Si ■ Su 


-0.01650 


0.12180 


-0.04887 


0.36408 


0.45061 


0.73273 


Si ■ Sis 


0.11611 


0.02574 


0.47814 


0.10860 


0.56423 


0.50655 




FIG. 1: Projection of the dodecahedron on a plane. The solid lines are antiferromagnetic bonds 
J. The black circles are spins Sj. 
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FIG. 3: Specific heat C of the antiferromagnetic Heisenberg model as a function of temperature 
T. Solid line: icosahedron with S{ = \, dotted line: dodecahedron with Sj = |, dashed line: 
icosahedron with Sj = 1. (7 is in arbitrary units and T in units of energy. 
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0,1 



1 

T 



10 



FIG. 4: Magnetic susceptibility \ °f the antiferromagnetic Heisenberg model as a function of 
temperature T. Solid line: icosahedron with Sj = ^, dotted line: dodecahedron with Si = ^, 
dashed line: icosahedron with s« = 1. x is in arbitrary units and T in units of energy. 
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1 



0.8 



0.6 - 



0.4 



0.2 



FIG. 5: Ground state magnetization M as a function of magnetic field h for the dodecahedron. M 
is the total spin S normalized to the number of sites and the magnitude of spin. Solid line: Sj = |, 
dotted line: Sj = 1. M has no units and /i is in units of energy. The steps between S sectors are 
0.1 for Si = \ and 0.05 for Sj = 1. The discontinuities are between 0.4 and 0.6 for s\ = ^, and 
between 0.4 and 0.5, and 0.7 and 0.8 for Sj = 1. 
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